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It is shown that low Reynolds number fluid flows can cause suspended particles to re- 
spond as though they were in an equilibrium system with an effective potential. This 
general result follows naturally from the fact that different methods of moving parti- 
cles in viscous fluids give rise to very different long-range flows. Two examples are dis- 
cussed: electrophoretic 'levitation' of a heavy charged sphere, for which a hydrodynamic 
'pseudo-potential' can be written in closed form, and quasi-two dimensional crystals of 
like-charged colloidal spheres which form near charged walls, whose apparent attraction 
arises not from a force but from persistent fluid flows. 



1. Introduction 

It is well-known that a particle which is moving under the influence of an external 
force such as gravity in a viscous fluid sets up a disturbance flow in the fluid which, in 
the low Reynolds number limit, decays with distance like r^^. The long-range character 
of these force-driven flows manifests itself in many situations. One example involves 
collective diffusive behavior. A pair of Brownian spheres has been predicted (Batchelor 
1976) and measured (Crocker 1997) to diffuse in a correlated manner. On the basis of 
these long-range correlations, Crocker et al. (2000) have developed a new technique for the 
microrheology of complex materials. Recently, Dufresne et al. (2000) experimentally and 
theoretically showed the diffusive behavior of a pair of Brownian spheres to be strongly 
influenced even by distant solid boundaries. Long-range flows also have a significant 
influence on the properties of sedimenting systems, giving rise to numerous divergent 
integrals (see e.g. Hinch 1988). For example, Caflisch & Luke (1985) have predicted 
velocity fluctuations in sedimenting particulate suspensions to diverge with system size, 
and Brenner (1999) has drawn attention to the fact that sedimenting particles' long-range 
interaction with cell walls can screen the flows and cut off the divergent fluctuations. 
In addition, long-range hydrodynamic coupling has provided a plausible interpretation 
(Squires & Brenner 2000) for the experiment of Larsen & Crier (1997), which measured an 
apparent attraction between like-charged colloidal spheres in the proximity of a similarly- 
charged wall. 

It is also well-known that there are other ways to move particles which do not give 
rise to such long-range disturbance velocity fields. One example is electrophoresis (e.g. 
Anderson 1989), wherein charged particles move under the influence of an applied electric 
field. A charged particle in solution attracts oppositely charged counter-ions from the 
solution, which form a screening 'double-layer' (e.g. Russel et al. 1989, chap. 4). The 
particle/double layer system is electrically neutral, so an applied electric field exerts no 
net force on the ensemble. However, the charged particle and the counter-ions are driven 
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in opposite directions, giving an apparent 'slip velocity' at the edge of the double layer, 
and a net motion to the particle. Since velocity fields with decay arise only from a 
net force on the fluid, electrophoretic (and all force-free) flow fields must decay at least 
as fast as r~^. 

In this paper we will demonstrate that this difference between the long-range distur- 
bance flows due to forced and force-free motions can be exploited to cause particles to 
behave as though they were subject to equilibrium forces derived from an effective po- 
tential. As an example, we will discuss the electrophoretic 'levitation' of a heavy charged 
sphere off a wall, for which an effective pseudo-potential 4>ps can be (lctc;rmined ex- 
actly. We will then present a method for making like-charged colloidal spheres near a 
wall behave as though there was an apparent attraction between them, so that layered 
like-charged colloidal crystals can bo grown. Throughout this work, we use the term 
'pseudo-potential' to emphasize that the hydrodynamic component of $ps arises not 
from a force, but from persistent low- Reynolds number flows which entrain each particle. 

We shall first briefly review the liydrodynaniics of forced- and force-free motion. Col- 
loidal particles (~ 1/xm) and their typical velocities l^m/s) are small enough that 
inertial effects are entirely negligible. In this limit, a rigid sphere of radius a moves due 
to an external force F at velocity Uo = boF, where 60 = (BTrrya)"^ is the Stokes mobility. 
This motion sets up a disturbance flow in the surrounding fluid 

"(^) = 8^U + l?rj-F-^(-^ + ^j-F, (1.1) 

where 77 is the fluid viscosity and R = r — is the vector connecting the observation 
point r and the center of the sphere tq. This can be re-expressed using fundamental 
singularities of the Stokes equation as 

The first term represents the flow due to a point force, known as a Stokeslet, which decays 
like R~^. The second term is the Green's function for a point source dipole, also called a 
potential dipole, which decays like R~^. The flow far from a forced sphere is dominated 
by the Stokeslet, which depends on neither the size nor the velocity of the sphere. 

By contrast, a charged insulating sphere with a thin double-layer moving due to a 
uniform electric field E^o sets up a disturbance flow outside the double layer (e.g. Russel 
et al. 1989, p. 256), 

eC f I 3RR\ ^ 

where e is the dielectric constant of the fluid and C is the potential difference across the 

double-layer. Here the sphere moves at velocity Vq = AfoEoo, where Mq = eC/47r?7 is 
the electrophoretic mobility. This flow is exclusively potential dipole flow and decays like 
R-^ 

Walls affect these flows significantly, since the fluid velocity must vanish identically on 
solid boundaries. Image singularities, analogous to image charges in electrostatics, can 
often be found which exactly cancel a flow set up on the wall. In this paper, we shall 
make extensive use of the flow due to a Stokeslet (point force) Fz oriented perpendicular 
to a wall, which Blake (1971) found to be 



u(R) = ( S(R) - S(R') + 2/1^ 
87r?7 



{v.v.(i)} 



- 2h 



dSjr) 

dzo 



(1.4) 
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where = R + 2hz is the vector between the images and the observation point, deriva- 
tives with respect to the source point ro (denoted by Vo) should be taken first, and 
then the position R' inserted in place of r. The flow due to the images decays hke R~^, 
whereas the image flow for a source dipole (Blake & Chwang 1974), which will not be 
required in this paper, decays like R~^. 

The velocity of a sphere forced toward a wall has two contributions: it moves at Uq due 
to the applied force, and it is advected with the flow set up by the image singularities. 
Evaluating the image flow from (1.4) at ro results in an 0{a/h) wall correction to the 
(scalar) mobility for motion perpendicular to the wall (e.g. Happel & Brenner 1983): 



h^{h)=bo 1 



9a ^ 
8^ + ^ 



/l3 



(1.5) 



The correction to the clcctrophorctic mobility can be derived in similar fashion, but is 
slightly more complicated because the wall gives both hydrodynamic and electrostatic 
contributions. Since the image flow (and electrostatic image potential) for electrophoretic 
motion decays like R^^ , the wall correction to the clcctrophorctic mobility is 0{a? /h^), 
and has been calculated using the method of reflections by Keh & Anderson (1985) to 
give 



8^3 



/l5 



(1.6) 



Throughout this paper, we will treat a/h and a/r as small, and so neglect this 0{a^/h^) 
variation in M±. 



2. Electrophoretic 'levitation' 

These two motions and their respective flows can be superposed, since Stokes' flow is 
linear. A charged sphere which sediments due to a force Fg towards a planar electrode 
and which is electrophoretically driven upwards due to a uniform electric field E^o will 
have a velocity v given by 

V = -b±{h)Fg + M^{h)Eoo, (2.1) 

with 6_L(ft) and M±{h) given by (1.5) and (1.6) respectively. If the bulk electrophoretic 
velocity is less than the bulk sedimentation velocity, so that ip = MoE^/boFg^ < 1, there 
will be a unique height 

where the two velocities exactly balance and a non-Brownian sphere will come to rest. A 
sphere started above ho falls faster than electrophoresis drives it up, and it moves down 
towards the plane. Similarly, a sphere below ho falls too slowly, and electrophoresis drives 
it up. Thus ho is a stable steady location. On the other hand, ho would be an unstable 
steady point for a sphere falling away from a wall. 

A sphere undergoing Brownian motion will have a probability distribution P{h) of 
being located at h which is given by the steady Fokker-Planck equation (see e.g. Batchelor 
1977) 

dP 

_ = = -V • J = -V • [v(/i)P - kBTh{h) ■ VP] , (2.3) 
where fcBT'bisatensordiffusivity. This equation has solution P(/i) = Pq exp(— ^ps/fcsT), 
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where $ps is a non-equilibrium 'pseudo-potential' 

^Mh)=F,h-f^^^dh'. (2.4) 

Inserting the values of b±_ and Mj_ from (1.5) and (1.6) and keeping terms to 0{a/h), we 
integrate (2.4) to obtain 



9a, fh 9 
■ In 

a 



(2.5) 



which can be expanded about its minimum at Hq to give an approximately harmonic well 



$ps 8 Fga 



(1 - v^)2 



jh h,f $ps(/^o) 

2a2 + ~fc^- ^^-^"^ 



If one viewed only the motion of a single sphere in this configuration, one would not 
be able to distinguish between this non-equilibrium pseudo-potential and a true ther- 
modynamic potential. However, the physical distinction is crucial: the only force on the 
sphere/double layer system is Fg, directed toward the wall. 

3. Multi-particle pseudo-potential interactions 

Even though the electropliorctically 'levitated' particle is stationary, it drives a long- 
range flow. Far from the sphere, the flow is always dominated by the force-driven motion, 
since it decays the most slowly. Streamlines for the flow around a hovering sphere, de- 
scribed in section 2, arc shown in figure 1. 

A second identical sphere, placed at the same height and at some distance r away from 
the original sphere, will both set up a long-range flow as well as be entrained by the 
flow set up by the first sphere. Forcing two spheres perpendicular to the wall induces a 
relative velocity between them (Squires & Brenner 2000), in the plane of the wall, as is 
evident from the streamlines in figure 1. This relative velocity of sphere 1 arises from its 
entrainment in the fiow set up by the images of sphere 2, and is found to 0{a/h) and 
0{a/r) from (1.4) to be 

rff_ 

27r?7 (r2 + 4/i2)5/2 

where x is the coordinate parallel to the wall, along the line joining the spheres, as in 
figure 2. Two spheres forced away from a wall move together as they leave the wall, 
whereas two spheres falling towards a wall move apart. In this case, we observe that 
the relative motion between two non-Brownian hovering spheres resembles the motion 
due to an effective repulsive force, even though it does not arise from a force. Example 
trajectories for two such spheres are shown by the solid curves in figure 2, with arrows 
showing the collapse onto these trajectories from various initial configurations. 



xi = -±2 = -— Fg, (3.1) 



4. Like-charge colloidal crystal formation 

Finally, we explore the possibility of setting up an attractive interaction between the 
spheres, which can lead to crystal formation. Since forcing two spheres toward a wall led 
to an effective repulsion, forcing them off of the wall will give an effective attraction. One 
way to do this would be to clcctrophorctically drive neutrally buoyant charged spheres 
into a similarly-charged wall. This would require the wall to possess a static surface 
charge density and to simultaneously passes a current. 
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Figure 1. Streamlines for a heavy charged sphere electrophorctically levitated at ho = 4a. Flows 
move down and radially out along the streamlines. The closed streamline shows the slip velocity 
on the sphere surface. The flow decays slowly and asymptotically approaches the flow due to a 
falling sphere near a wall, even though the sphere is immobile. Gravity exerts a force directed 
toward the wall, and given a sphere with positive C potential, the electric field is directed away 
from the wall. 



As wc have seen, the far-field hydrodynamic flows are dominated by the net force 
on the body and its charge cloud. In this case, the net force F^^ on each is given by 
its electrostatic repulsion from the static charge on the wall, found by integrating an 
electrostatic stress tensor over the surface of the sphere (see, e.g. Russcl et al. 1989). 
Under the linear superposition approximation, the repulsive force the wall exerts on the 
sphere is found to be 

/^=47^ZagAB-^e-«^ (4.1) 

where Z is the effective charge on the sphere, Ug is the charge density on the wall, 
is the Debye screening length and Ab is the Bjerrum length(?a .7 nm in water at room 
temperature). This force sets up a flow around each sphere like in figure 1, with flow 
arrows reversed. The lateral component of the flow now draws the spheres together. 

Two non-Brownian spheres separated by distance r are immobilized at the height ho 
where M • = h{ho,r)-Fy,{ho). The equivalent force on each sphere which would be 
necessary to drive this motion in a quiescent fluid is 

« , 1\ . (4-2) 

where 6x1x2 is the coefficient of b which gives the mobility sphere 1 in the x direction 
due to a force on sphere 1 in the same direction, and so on. Using (3.1), and neglecting 
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Figure 2. Two heavy spheres electrophoretically driven upwards approach a stable trajectory, on 
which the spheres move apart. This mimics a repulsion between the spheres, and asymptotes to 
two spheres levitating independently. In this figure, we have chosen E such that an independent 
sphere would hover at ft = 6a. 



the weak dependence of ho upon r, this can be integrated to yield a pseudo-potential 

r(4/i2+r2)-5/2 



$ps(r, ho) « Up{r) 



5hlF^{ho) 



-dr. 



(4.3) 



27r77 J bxixiiho) -bj:^a:2{rJ'-Q) 

where Up is the true pair potential between the spheres. Although. bx2Xi 

is 0{a), in 

practice it is small compared to fexixn and its neglect allows the direct integration of 
(4.3) to give 

The true pair potential between charged spheres, Up, is given in linear superposition 
approximation by (e.g. Larsen & Grier 1997), 



U 



" =Z^Xb 



1 + Ka 



(4.5) 



Equation (4.4) is identical to that derived by Squires & Brenner (2000) in their analy- 
sis of the experiment of Larsen & Grier (1997) which measured the effective interaction 
between like-charged colloidal spheres near a single charged wall. Although the two for- 
mulae were derived in the same fashion, they have very different physical consequences. 
In that experiment, a pair of spheres was repeatedly released from the same height ho and 
then re-trapped; therefore, a net force systematically drove the spheres away from the 
wall throughout the experiment. Consequently, there was a systematic relative motion 
of the spheres towards each other which had been interpreted to arise from an attractive 
force. In an equilbrium system, where the spheres fluctuate into the wall as often as away, 
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one would expect no systematic relative motion, and this hydrodynamic effect would be 
of no physical consequence. 

Similarly, the pseudo-potential of equation (4.4) does not describe an attractive force 
between two like-charged colloids. However, in the present situation, similarly- charged 
spheres actually do respond as though there was an attraction. The effect of electrophore- 
sis is to 'pin' the spheres at some height without applying a force, so that the long-range 
flow is dominated by the forced motion and entrains neighboring spheres accordingly. A 
collection of spheres can thus aggregate to form clusters-a behavior which might other- 
wise seem to arise from an attractive force. 

Spheres which are not neutrally buoyant have a similar pseudo-potential, with Fy, — Fg 
in place of Fw in (4.4). Note that both the depth and the range of this potential can 
be tuned independently-the range by adjusting the electric field Eoc, to give different 
heights ho and the overall magnitude of the attraction by adjusting the relative density 
of the particle to affect Fg. Furthermore, one could study the freezing and melting of the 
crystals we'll describe next by simply adjusting Eoo during an experiment. 

Brownian motion considerably complicates the above picture. Again, the probability 
distribution for the spheres' positions obey the steady Fokker-Planck equation (2.3). An 
interesting fact is that there is a persistent non-zero probability current J in the steady 
state. The probability current can only be zero when 

-V$ps = • u (4.6) 

has a solution. In equilibrium systems, u = b • F so that F = — V$ps as expected. Here 
u = b • F -|- u_E, and (4.6) will not in general have a solution. Physically, this is because 
the equivalent force required to give velocity ue would be non-conservative. Vertical fiuc- 
tuations take the pair into regions of stronger or weaker attraction, leading to persistent 
closed orbits of J. When these vertical fluctuations are relatively small, J is small and the 
effective pseudo-potential is well approximated by (4.4). A comparison betweeen (4.4) 
and simulated results is shown in figure 3. The lower two pseudo-potential wells (a) and 
(b) are deep enough that bound states form. The simulated pseudo-potentials ^pg/ksT 
are thus simply — log g{r), where g{r) is the the pair correlation function obtained from 
simulations of a pair of spheres. We used the method described by Crocker & Grier (1996) 
to simulate the shallowest potential well (c), since it is not deep enough for bound states. 

We have run simulations which account for Brownian motion to demonstrate that 
spheres can indeed self-assemble into a quasi-2D crystal. We use the dynamical equations 

v = b • F + M • E (4.7) 

where v,b and so on are multiparticle quantities. We retain terms up to 0{a) in b and 
M and use a box of side lOOo with periodic boundaries. Brownian motion is accounted 
for in standard fashion (Ermak & McCammon 1978: Grassia et al. 1995), wherein D = 
ksTh is a position-dependent tensor diffusivity. We use as input parameters the values 
Z = 7300, = .275^rn and 2a = .652^rn as in the experiments of Larsen & Grier 
(1997) and the value o"g = Adg for the (unnic;asured) charge density on the wall which 
Squires & Brenner (2000) found as the best fit to the experimental data. We choose 
these values because 1) the effects we are discussing are most pronounced with long 
screening lengths and high surface charge densities, so that the spheres tend to 'hover' 
at an appreciable distance from the wall, and 2) these effects may bear some relation 
to these particular experiments. We started nineteen spheres in a hexagonal crystal, 
then simulated their dynamics until it was clear that the qualitative features would not 
change. This was typically approximately a minute or less of real time. Figure 3 shows 
the resulting suspensions at the end of each run for three different field strengths. The 
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strongest field, for which ho = 2/xm, yields a hexagonal crystal, whereas a weaker field, 
for which ho = 2.2/zm, is strong enough to yield a dense phase, but not strong enough 
to retain orientational order. We loosely term this an istropic 'liquid' phase. Finally, for 
an even weaker field, for which ho = 2.5/im, the attraction is sufficiently weak that the 
spheres spread out to fill the box as an isotropic 'gas'. 




2 3 4 5 6 r(|i) 7 

Figure 3. Radial pseudo-potentials and representative 'phases' of neutrally buoyant charged 
spheres electrophoretically driven into the double-layer surrounding a charged wall, a) A strong 
field (ho = 2/im) yields a deep pseudo-potential well and tight pairwise confinement, giving a 
regular triangular lattice, b) A moderate field {ho = 2.2fj,m) yields a pseudo-potential well deep 
enough to yield a dense phase, but not strong enough to retain orientational order, c) A weak 
field [ho = 2.5pm) yields an isotropic 'gas'. Dotted lines represent 1 ksT of energy and give 
typical excursions about the equilibrium point. Lines give the approximate analytic form (4.4), 
and dots give the simulated pseudo-potentials. 



5. Discussion and conclusions 

While wo have specifically investigated the combination of clcctrophorctically-drivcn 
motion and body-force motion, the present results have a wider validity. Our results 
will hold whenever motion using body forces such as gravity, electrostatic repulsion or 
attraction, magnetic forces, etc. arc combined with force-free motion. There arc many 
ways to achieve the latter. Since any self-propelled motion can not result in a net force 
or torque on the system, all swimmers generate flows decaying faster than (Lighthill 
1976). There arc also other types of "phorctic" motion, which involve particle motion 
due to fields which interact with the particle's surface (e.g. Anderson 1989), such as 
thermophoresis (motion due to thermal gradients) and diffusophoresis (motion due to 
solute concentration gradients) . Similarly, the thcrmocapillary motion of bubbles (Young 
et al. 1959), driven by thermal surface tension gradients, is force- free. 

Having shown that like-charged colloidal crystals can self-assemble into ordered crystals 
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due to long-range persistent non-equilibrium hydrodynamic flows, it is natural to ask 
whether there could be any connection between these crystals and those observed by 
Larsen & Grier (1997). These metastable crystallites formed near the charged glass walls 
of the experimental cells, and persisted for times orders of magnitude longer than purely 
repulsive Brownian spheres should. This response can certainly be understood within the 
present picture, and would require a force- free motion to drive the spheres into the charge 
cloud of the walls. A spurious current in the cell could do this; however, it is unlikely 
that this explains Larsen and Grier's crystals, for two reasons: 1) there are no electrodes 
on the top and bottom walls of their cells, and 2) the crystals were observed to form 
on both the top and bottom walls, whereas a net current in the cell would likely drive 
the spheres to one wall or the other. Another possible current source, suggested by D. S. 
Fisher, is the dissociation of ionic groups on the glass surface. Other possible mechanisms 
for force- free motion include thermal or electrolyte concentration gradients; however, it 
is not clear which if any of these mechanisms are present or responsible for the crystals. 
It is generally believed that a novel long-ranged attraction between confined like-charged 
colloids, whose origin remains controversial, gives rise to this behavior (e.g. Hansen & 
Lowen 2000). The present work provides a non-equilibrium mechanism whose origins are 
well-known and understood which would mimic an attractive interaction and could drive 
like-charge colloidal crystallization, so its potential role should be investigated further. 

Two-dimensional colloidal crystals have been formed on electrodes via electrophoretic 
deposition (Bohmer 1996; Trau et al. 1997), and Solomentscv et al. (1997) has proposed 
a electrohydrodynamic model, which bears some similarity to the present work, for its 
explanation. The crucial difference between the our work and that of Solomentsev et al. 
(1997) and is that in the latter, the particles arc first deposited on the electrode, and then 
electro-osmotic flows drive adjacent particles towards one another. These electro-osmotic 
flows are much shorter-ranged than the body-force driven flows described here, and are 
highly screened by the close proximity of the wall. 

In conclusion, we have used the difference in flow fields set up by force-driven mo- 
tion and force-free motion to show that persistent long-range viscous flowflelds can be 
set up, resulting in particle motion which mimics an effective potential. This appears 
surprising at flrst, since the system is out of equilibrium, and yet behaves as though it 
is an equilibrium system with an effective potential which is hydrodynamically driven. 
Given the general nature of this effect, these persistent flowflelds are likely to flnd various 
applications in self-assembly and microfluidic contexts. 

This work was supported under the NSF Division of Mathematical Sciences grant 
DMS-9733030. 
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